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AMENDMENTS TO THE SPECIFICATION : 

Before the paragraph beginning at page 1, line 4, 
insert the following : 

— REFERENCE TO COMPUTER PROGRAM LISTINGS 

A computer program listing appendix with twelve 
files is provided on a compact disc as part of the invention 
disclosure. The material of the twelve files on the compact 
disc is incorporated by reference. The files included on the 
compact disc are: 

• generate-discretization.mu has a size of 6,241 bytes and was 
created (stored on the CD-R) on December 12, 2006. 

This file contains the program "generate-discretization" , 
referred to as Appendix 1. 

• preparations .mu has a size of 3,101 bytes and was created 
(stored on the CD-R) on December 12, 2006. This file contains 
the subroutine "preparations", referred to as Appendix 2. 

• setup-equations .mu has a size of 17,414 bytes and was 
created (stored on the CD-R) on December 12, 2006. 

This file contains the subroutine "setup-equations", referred 
to as Appendix 3. 

• solve-equations .mu has a size of 4,681 bytes and was created 
(stored on the CD-R) on December 12, 2006. This file contains 
the subroutine "solve-equations", referred to as Appendix 4. 

• analyze-solution.mu has a size of 22,654 bytes and was 
created (stored on the CD-R) on December 12, 2006. This file 
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contains the subroutine "analyze-solution" , referred to as 
Appendix 5 . 

. appendix6.txt has a size of 10,588 bytes and was created 
(stored on the CD-R) on December 12, 2006. This file contains 
the output for D100, Order 2, on the grid 
(-1,1) "3, optimize=0, referred to as Appendix 6. 
. appendix7.txt has a size of 3,289 bytes and was created 
(stored on the CD-R) on December 12, 2006. This file contains 
the output for D100, Order 2, on the grid 
(-1,1) "•3, optimize=l, referred to as Appendix 7. 
. appendix8.txt has a size of 9,465 bytes and was created 
(stored on the CD-R) on December 12, 2006. This file contains 
the output for D100, Order 2, on the grid 

(-1,1) A 3, optimize=2, referred to in the following as Appendix 
8. 

• appendix9.txt has a size of 47,206 bytes and was created 
(stored on the CD-R) on December 12, 2006. This file contains 
the output for D100, Order 4, on the grid 
(-2,2)^3, optimize=0, 30, referred to as Appendix 9. 
. appendixl0.txt has a size of 4,092 bytes and was created 
(stored on the CD-R) on December 12, 2006. This file contains 
the output for D100, Order 4, on the grid 
(-2,2) A 3, optimize=l, referred to as Appendix 10. 
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• appendixll.txt has a size of 49,617 bytes and was created 
(stored on the CD-R) on December 12, 2006. This file contains 
the output for D100, Order 4, on the grid 

(-2,2)^3, optimize=2, 5, referred to as Appendix 11. 

• appendixl2.txt has a size of 16,572 bytes and was created 
(stored on the CD-R) on December 12, 2006. This file contains 
the output for D200, Order 2, on the grid 

(-1,1) A 3, optimize=0, referred to as Appendix 12.-- 

Please replace the paragraph spanning pages 8 and 9 
with the following: 

--An example is a general one-dimensional discretization for a derivative at node 
i which uses the stencil between the nodes i - m and i + n, where m and n are given positive 
integers. On the stencil S = (i - m,i -m + - l,x,i + + n - 1, i + n) , the approximation 

of the first derivative can be written as 

1 r 

u x =-—{a^ m u i _ m + a_ m+1 Ui„ m+1 + ... 
A* 

+ fl-lMi-l + W + a l U i+l + * • 
a n-l u i+n-l + a n u i+nl 

0) 

The coefficients fly , i - j e S fly , j e S determine the nximerical properties of the discretization. 

A general description of the above discretization can be found in Dl which discusses the 
method for obtaining the coefficients fly. Expression (1) also generalizes to higher 

derivatives. — 

Replace the paragraph beginning at page 47, line 16 
with the following: 
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2 2 

-The second derivative D 2 - 9 * D 2 = can be expressed in the grid-based 
derivatives u xx ,1*^, ,u z2 ,w xy ,u y2 and tt 2X according to 



t» = cos 2 a cos 2 /^ + sin 2 a cos 2 >^ + sin 2 /*^ + 
del dy dz (n) 

2cosacos 2 >3sina|^ + 2cosacos/?sin/?|^ + 2sinacos > 0sin/?^ 

The terms T xx ,T xy ,T xz ,T yy ,T yz and T 2Z are obtained by the computer program mentioned 
before, and are given in appendix 12. The stencils T xx ,T xy ,T XZ ,T yy ,T yz and T zz have been 
added in various quantities to the stencils Tf^TfoTf^TfoTfo and T/ 2Z to obtain the 
more symmetric representation of equation 20. This shows once more the use of the degrees of 
freedom, and the equivalence between two expressions for an approximation of D p using a 
different Tf but sharing the basis described by the stencils.- 

Replace the list of symbols beginning on page 54 with 
the following: 

—List of symbols 

A an amplitude of a Fourier component 
A', A" intermediate bases in the transformation from the grid basis 
to the local basis B 
a,b,c,. . . components of the vector a 
a vector of preferential direction 
B(h ,e 2 ,e 3 ,..)a local basis, with e x along a preferential direction, i.e. e a / / a 
C c computational coefficients used in the approximation of D p 



which are dependent on the numerical formulation used 
C s computational coefficients used in the approximation of D p , 

in the Finite Difference formulation ; weighting coefficients 
"/ m n wei gh tin g coefficients C s for node /, m, n, • • 

D p spatial p th derivative, to be discretized 
Dp an approximation to D p 

D p C an approximation to D p in the Finite Difference formulation, 
representing a linear combination of values 
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Dp 1 an approximation to D p in the Distribution Method, 
depending on the distribution coefficients a z 
D^ f¥ an approximation to D p in the Finite Element formulation 

depending on the basis function <p and the test function i// 
f the flux 

2 

J the imaginary unit, such that I = — 1 

if j, k,-' indices numbering the nodes of a structured grid 
iax / /max ' *max maximum indices of a grid 

I e l the integral of the derivative D p over a volume, used in the 

Residual Distribution Method 
M order of the error of a discretization 
N number of dimensions 
P the point where the derivative is computed 
p index for a higher ( p th ) derivative, or first derivative ( p = 1 ) 
Pi ' P2 f P3 * ' " ^ e P owers of the derivatives with respect to , e 2 , . . • in a 

mixed derivative 
qi , q>i , ^3 , • • • arbitrary variables 

r an integer summation index 
r max th e maximum value of r in the summation 

S the stencil : the set of points used in the computation 
of the approximation D p 

t the time coordinate 

' 11 ' *1 2 / " coefficients used in the transformation between coordinate 

systems 

T p represent the terms in the discretization f3 resulting from 

degrees of freedom which remain when the approximated value 
Dp is optimized 

u unknown at a grid point 
u s unknown at a point of the stencil S 

du 

u a derivative of u with respect to a , e.g. u a = — — 

doc 

x N -dimensional position vector 
V differential operator, , ^ , , • • J 
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a { the distribution coefficient used in the Distribution Method 
for the distribution of the part a { l ei or a\f to node i 
Ax, Ay, Az, • • the mesh spacings in the coordinate directions 
At the time increment 

— partial p th derivative wife of u with respect to x 

dx 
gP 

derivative with respect to x 



dx? 

8 n error term in D p 

s s error term in the expression of u s using a truncated 

Taylor series expansion 
k the wave number vector 

tp the basis function used in the Finite Element method for repre- 
senting u over the element 
y/ test function used in the integrals of the derivative in the Finite 

Element formulation 

co the angular frequency- 
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